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Abstract 

We examine effects of encapsulated metal atoms inside a Cqq molecule on the activation energy 
barrier to the Stone-Wales transformation using ab initio calculations. The encapsulated metal 
atoms we study are K, Ca and La which nominally donate one, two and three electrons to the 
Ceo cage, respectively. We hnd that isomerization of the endohedral metallofullerene via the 
Stone- Wales transformation can occur more easily than that of the empty fullerene owing to the 
charge transfer. When K, Ca and La atoms are encapsulated inside the fullerene, the activation 
energy barriers are lowered by 0.30, 0.55 and 0.80 eV, respectively compared with that of empty 
Cqq (7.16 eV). The lower activation energy barrier of the Stone- Wales transformation implies the 
higher probability of isomerization and coalescence of metallofullerenes, which require a series of 
Stone- Wales transformations. 

PACS numbers: 81.05.Tp, 34.10,+x, 61.48. +c 
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Since its discovery, the fullerene has been extensively studied both theoretically and ex- 
perimentally because of its potential applications in the fields of nanomaterials and biomed- 
ical science. Although the fullerene molecule usually has the highest icosahedral symmetry 
allowed by sixty constituent carbon atoms, transformations to other isomers of lower sym- 
metry are still possible. One of the most plausible processes for isomerization of the fullerene 
is the so-called Stone- Wales (SW) or "pyracylene" rearrangement, which is the 90° rotation 
of two carbon atoms about the midpoint of the bond. Q] The SW transformation is also used 
to describe the structural changes of sp 2 -bonded carbon nanosystems. JjJ For example, it has 
been proposed that the fusion process of fullerenes or carbon nanotubes may be explained 



in terms of a sequence of this rearrangement.^, |5j,|6( A problem in this proposal is that the 
activation energy barrier of ~ 7 eV,|7J on the basis of cutting two a bonds, is somewhat too 
high to overcome at experimental temperature for the merging process (1000— 1500°C). 
Nevertheless the transformation is useful to explain the production of all isomers of fullerenes 
and carbon nanotubes. Eggen et al. found that the activation barrier for this rearrangement 
was substantially reduced in the presence of a loosely bound carbon atom located prefer- 
entially in the region of paired pentagons. [loT] Therefore pentagon rearrangements which 
are necessary steps in the growth of fullerenes may considerably result from autocatalysis 
by carbon. Osawa et al. reported the catalytic effect of various kinds of elements or ions 
in the bowl-shaped C^Hn.|llj However, relatively few studies have been devoted to the 
SW transformation of the endohedral metallofullerene. In this paper, we investigate effects 
of encapsulated metal atoms on the activation energy barrier in the SW transformation of 
the M@C6o molecule, where M = K, Ca or La. We shall show below that isomerization of 
the endohedral metallofullerene via the Stone- Wales transformation can occur more easily 
than that of the bare fullerene owing to the charge transfer from the encapsulated metal 
atoms to the fullerene. Experimentally, metallofullerenes are much easier to fuse than empty 
fullerenes in peapods. |3, |l^ 

In the present study, we perform ab initio pseudopotential calculations using the plane 
wave basis setjlJ] within the local den sity approximation (LDA) for the exchange-correlation 
effect. The ultrasoft pseudopotential|15| is adopted with the cutoff energy of 30 Ry. The 
supercell size is 29x29x29 a 3 B , where as is the Bohr radius, which can keep the undesir- 
able interfullerene interaction sufficiently small. For La, we treat all 11 valence electrons 
(5s 2 5p 6 6s 2 5(i 1 ) explicitly since the 5s and 5p orbitals cannot be considered to be frozen as 
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core electrons. Q] The scalar relativistic effects are taken into account because the effects 
are important in heavy elements such as La. We start the optimization of the La position on 
a three-fold axis 1.5 A away from the center of the cage and obtain the stable position of the 
La atom in the Cm cage. For K@Ceo and Ca@Ceo, the starting positions for optimization 



are at the center 



17] and on a five- fold axis 1.5 A away from the center of the cage,|l8fl 



respectively. The atomic positions are relaxed until the normal component of the atomic 
force becomes less than 0.02 eV/A. We use the nudged elastic band (NEB) method[l^ to 
calculate the activation energy barrier in the SW transformation of the fullerene. The NEB 
method is efficient in finding the reaction pathway with minimum energy when both the ini- 
tial and final states are known. Regarding SW transformation of metallofullerene, we rotate 
the C-C dimer which is the closest to the metal atom (Ca and La). Thirteen replicas are 
chosen including the initial and final ones to construct the elastic band. We also calculate 
the climbing image of the highest energy replica to pick up the exact transition state which 
is regarded as a saddle point in the potential energy surface. 

In the relaxation of the metal atom in Ceo, the K atom prefers to be at the center of Ceo- 
On the other hand, the stable position of the La atom is off-center and the distance between 
the La atom and the nearest C atom is ~ 2.5 A. The Ca atom shows the same tendency 
as La. Upon encapsulation of the La atom, the three-fold degeneracy of the ti u states is 
split into 2+1 since the off-centered La atom breaks the icosahedral symmetry of an isolated 
Ceo- 20] Almost three electrons are transferred from the La atom to the doubly degenerate 

nn 

states (four- fold degeneracy including spins) of Ceo-|20,|21| The two-fold degeneracy is split 
again by the charging of odd number of electrons. It is found that a nearly symmetric 
pathway which has a C2 symmetry appears at the saddle point (transition state) and the 
activation energy profile [shown in Figs. 1(a) and (b)]. Comparing Ceo and La@Ceo, the 
rotating carbon dimer has a bond length of 1.23 (1.27) A, with 1.40 (1.39) A bonds to each 
of the two nearest carbon atoms for Ceo (La@Ceo)- This result is in agreement with the 



previous studies. |7J, 



17. 
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We identify the atomic geometry of the transition state (TS) in SW transformations of 
Ceo and La@Ceo as depicted in Figs. 3(a) and (b). Table I lists the activation barriers (E a ) of 
the SW transformations in various sp 2 -bonded carbon systems. E a in the SW transformation 
of the fullerene is lower than in the graphene layer (9.2 eV, the highest among the systems 
we study). Lowering of the barrier in fullerenes is due to the curvature effect associated 
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with the deviation from the sp 2 bonding. |6] When K, Ca and La atoms are encapsulated 
inside the fullerene, the activation energy barriers are lowered by 0.30, 0.55 and 0.80 eV, 
respectively compared with that of bare Cgo- Reduction of the activation energy is smaller 
than in the case of the autocatalysis and other covalent materials. jloL 11 1 To examine the 
major effect of the encapsulated metal atom on the lowering of the activation barrier, we 
perform calculations for fullerenes charged with one or two electrons. Note that we do not 
show the data for Cg^. The third donated electron is not bound to the fullerene. When 
a fullerene is charged, the activation barrier becomes lower than that of the neutral cage. 
The barrier reduction is approximately proportional to the number of donated electrons. 
E a in the SW transformation of Cg (Cgg) is almost the same as that of K@C 60 (Ca@C 60 ). 
Therefore, the main role of the incorporated metal atom is the electron donation. 

We find the reaction energy (E r = E pro d U ct — E reactant ) shows the same trend as E a in 
Table I. The C 2l ,-symmetry C 60 with two pairs of fused pentagons (obtained by the SW 
transformation) is 1.59 eV higher in energy than the K-symmetry Cgo- This value is in 
agreement with other group's results in the literature. [23] By contrast, for La@Ceo, E r 
is 0.3 eV, which is much smaller than that of Ceo- It represents that an isomer of the 
metallofullerene which has adjacent pentagons can exist with higher probability than bare 
Ceo- In the case of K@C6o, E r has an intermediate value. The K atom at the center of the 
cage does not move during the process of the SW transformation. By contrast, the La atom 
in the fullerene moves by 1.0 A toward the paired pentagons in the process of dimer rotation. 
The displacement of the La atom induces a large elongation of the Ceo cage compared with 
the empty Ceo molecule. The geometry of an empty fullerene is governed by the isolated- 
pentagon rule (IPR),0, 13 which says that the most stable structure is that in which all 
pentagons are surrounded by five hexagons. The IPR-violating cages such as Sc2@C66 and 
Sc 3 N@C 68 , however, are stabilized by the electron transfer between the encapsulated metal 
atoms and the carbon cage, which significantly decreases the strain energy caused by paired 
pentagons and thus stabilizes the fullerene. [3, 13 Owing to the 108° bond angles of the 
pentagon, the apex atoms in the fused (paired) pentagon regions of the empty fullerene have 
a hybridization very close to sp 3 and exhibit a dangling-bond-like character which cannot 
be hybridized with any other electronic state. [28] For metallofullerenes, on the other hand, 
the electron donation from the metal atom weakens the dangling-bond- like character. 

To clarify the role of the encapsulated metal atom, we have investigated the electronic 
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structures of C 60 and La@C 60 in detail. The projected density of states (PDOS) to the 6s 
and 5d orbitals of La for each reactant, transition state and product are drawn (not shown). 
It is observed that the energy levels of the La atom lie mostly above the Fermi level. It 
means that almost three electrons move from La to the fullerene in La@Ceo- In Fig. 2(a), 
the highest occupied molecular orbital (HOMO) at —0.2 eV corresponds to the dangling- 
bond-originated state as illustrated in Fig. 2(b). This state is above the HOMO of Ceo 
(i.e., the five-fold states originating from the HOMO of empty C 60 which are split-off now). 
Energy level which has the dangling bond character in La@C 60 is at —0.7 eV (HOMO— 2) as 
shown in Fig. 2(c) by arrow and the wave function in Fig. 2(d). Here we can see the three 
electrons coming from the La atom occupy the upper one and half level which is the lowest 
unoccupied molecular orbital (LUMO) and LUMO+1 state in the transition state of C 60 . 
The overall shape of two wave functions [Figs. 3(b) and (d)] is the same with small change 
originating from the perturbation of the La atom. Figures 3(c) and (d) show the PDOS into 
2j?-orbitals of two carbon atoms [labeled A, B and A', B' in Figs. 3(a) and (b)] which have 
the dangling bonds in the transition state. In the case of empty C 60 , the PDOS of A and 
B are exactly the same. For La@Ceo, on the other hand, the charge density of nonbonding 
states at site A' (between —1 and eV) is not identical to the density at site B'. This 
difference arises from the hybridization between the La atom and the C atom at site B. In 
the TS configuration of empty C 60 , the HOMO and LUMO are both localized states within 
the hemisphere containing two dangling atoms (A and B). By contrast, the TS geometry 
of La@C 60 has the HOMO and LUMO delocalized on the whole cage. Consequently, the 
electron donation and hybridization between the fullerene and metal atom inside reduce the 
activation energy barrier of the SW transformation in the endohedral metallofullerene. 

We have studied the effect of encapsulated metal atoms inside a Cqo molecule on the 
activation energy barrier of the SW transformation using the ab initio density functional 
calculations. The metal atoms in our study are K, Ca and La, which donate one, two 
and three electrons to the C 6 o cage, respectively. We have found that isomerization of the 
endohedral metallofullerene by the SW transformation can occur more easily than that of 
the empty fullerene owing to the charge transfer. The reduction of the activation energy 
barrier by the electron transfer may also explain the fact that metallofullerenes are much 
easier to coalesce than empty fullerenes inside nanotubes. 
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[Table Caption 

Table I. Activation barrier (E a ) and reaction energy (E r ) of the SW transformation in 
various sj? 2 -bonded carbon systems. All energy values are in eV. 
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[Figure Captions 



Figure 1: Schematic energy diagram and the atomic configuration along the reaction 
coordinate in the Stone- Wales transformation of empty C 60 . (a) Energy profile along the 
reaction coordinate. The dots represent the energies of the replicas in the SW transformation 
in the NEB method, (b) Atomic rearrangement from the lh isomer to the C2 V isomer of the 
empty Ceo molecule. 

Figure 2: (color online) Density of states (DOS) for C 60 and La@C 60 in the transition 
state [(a) and (c) respectively] and the isodensity surface plots of particular wave functions 
for C 60 and La@C 60 [(b) and (d)]. The energy levels of the wave functions in (b) and (d) 
are indicated with an arrow in (a) and (c), respectively. The values for the red and blue 
isodensity surfaces are ±0.0014 e/A 3 . 

Figure 3: (color online) Atomic structure of the transition state configuration for (a) C 60 
and (b) La@C 60 , and PDOS of two atoms [labeled A, B and A', B' in (a) and (b)] which 
have the remnant dangling bonds in the transition state for (c) C 60 and (d) La@C 60 . The 
defect states due to the dangling bonds and distorted bonding configurations are indicated 
by arrows. The red solid and blue doted lines are the PDOS decomposed into 2p-orbitals of 
A, A' and B, B' atoms, respectively. 
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System E a E r 

C 60 7.16 1.59 

K@C 60 6.85 1.07 

Ca@C 60 6.61 0.59 

La@C 60 6.36 0.30 

C 60 6.87 1.12 

C 2 m 6.59 0.67 

C 42 Hi 6 (for graphene) 9.20 3.07 

Table I 

W. I. Choi, G. Kim, S. Han and J. Ihm 
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